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We present a detailed study of the dynamic response of a ring of N equally spaced 
Josephson junctions to a time-periodic external flux, including screening current 
effects. The dynamics are described by the resistively shunted Josephson junction 
model, appropriate for proximity effect junctions, and we include Faraday's law for 
the flux. We find that the time-averaged I—V characteristics show novel subharmonic 
giant Shapiro voltage resonances, which strongly depend on having phase slips or not, 
on N, on the inductance and on the external drive frequency. We include an estimate 
of the possible experimental parameters needed to observe these quantized voltage 
spikes. 
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Motivated by significant advances in microfabrication techniques of Josephson junction 
(J J) arrays, many contemporary studies have focused on understanding the dynamic re- 
sponse of two-dimensional J J arrays [1-7]. One type of array with well-defined uniform 
properties is made of superconducting-normal-superconducting (SNS) proximity effect junc- 
tions. Several novel collective locked-in states have been discovered which are manifested in 
the current versus voltage (IvsV) characteristics as plateaus. These giant Shapiro steps are 
a macroscopic manifestation of diverse underlying collective vortex oscillating states [1-5]. 
Most recent studies, experimental and theoretical [1-5,7], have considered current driven 
two-dimensional J J arrays. In this paper we study the dynamic response of a ring of N 
equally spaced J J, referred to as a Josephson Junction Necklace (JJN), driven by a perpen- 
dicular external time-dependent flux <3> e (t), including self-induced magnetic field effects. The 
motivation of this study is two- fold; first, to understand the dynamic response of a J J array 
to a $ e (t), and second, to assess the importance of geometry in this response. As we show 
in this paper, the physics of the JJN is, in many respects qualitatively and quantitatively 
different from that of a current driven planar array. 

The JJN is shown schematically in the inset of Fig. I. This model is a natural extension 
of the standard rf SQUID to N junctions instead of one. The external time-dependent 
flux considered here produces a Faraday voltage V e (t) of the form, — ^$ e (£) = V e (t) = 
Vdc + cos 2ni>t. We also include the self-induced magnetic flux $j due to the screening 
currents, of importance in a proximity effect J J. We model the dynamics of the order- 
parameter phase of the JJN in terms of the resistively shunted Josephson junction model 
(RSJ) defined by the current 

/(*) = / c sin[0i - - V> M+ i] + 2^^[& ~ ~ ^v+il- (!) 

In this equation A M = (pi — <t>i+\ is the phase difference for the /x-th junction denoted as x 
in the inset of Fig. 1. I c is its critical current and R its shunt resistance, both assumed to be 
the same for all junctions in the necklace. Here ibi i+i(t) — — A-d£ — — — = — fit), 
where <3>(t) is the total flux through the ring of radius r, f(t) = ( < I ) e (t) + $ i (t))/ < I ) o — fe + fi, 
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and $o = — is the quantum of flux. The induced flux is given by $j = LI(t), where L 

is the geometrical inductance of the ring. It is convenient to write Eq(l) in dimensionless 

units with i(r) = I(t)/I c , r = t(2eRI c /h), v = V/(RI C ), £ = L^^, and v+ = " . In 

$o 2eRI c 

these units Eq.(l) becomes 

i(r) = sin[^([/(r)] fc - /(r))] + ~(W)] fe - /(r)). (2) 

In writing this equation we have used the fact that the current is conserved in the ring 
which, together with the periodic boundary conditions (p.b.c.) 4>n+i = <fii, implies that the 
phase difference A(p = ^[f(r)]k, where [/]*. stands for the nearest integer function defined as 
[f( T )]k = hit(/ + -). The specific form of the applied external flux considered in this paper 
is f e = $e/$ = ~(-fe + W2vrr + j^ffa sin 2-kv^t) , with /° = / e (r = 0). Substituting 
this expression for f e in Eq(2) we get 

^i(r) + |i(r) - sin[^[/(r)], + ^(-2*/° + v dc r + ^- sin2vr^r - «(r))]| 

= -j^( V dc + Vac COS + ^^l[/( r )]fc- ( 3 ) 

Eq(3) is the general equation describing the response of the JJN to a time-dependent external 
flux, and understanding the structure of its time-dependent solutions is the central goal of 
this paper. To find i{r) from Eq(3) we need to know [f(r)]k- As shown below, we can 
use the ground state properties (f e = /°) of the JJN to obtain solutions for [f(r)] k in 
the relevant physical regimes. We start then by considering the ground state properties 
of the JJN including inductive effects, which appear not to have been discussed before for 
iV > 2 ||. From the p.b.c. and current conservation it follows that the phase difference 
for each and all junctions is A M = (2irk)/N, with k = 0, 1, 2, 3, .., N — 1. When the self- 
inductance is neglected, the ground state energy per junction E g /N, in normalized units, 
* s €n = io^N = ~ cos lw (ife)k ~ fe)], with the corresponding normalized ground state 
current ijv obtained from zjv = — f^if = sm [^f([/e]fc — fe)]- For a periodic JJN we note 
the important symmetry of the ground state energy /° — > /° + 1 and f® — > /° — 1, which 
allows the analysis to be restricted to /°e[0, 1]. In Fig. 1 we show €n(£ = 0) (solid line) 



for iV = 10. The difference between curves / and II is that in II there is a discontinuity 
in e^v at f® = 1/2, whereas curve I is continuous. Physically, as /° increases from zero to 
/,? = (1/2) - the JJN can either follow curve I or curve II. Following II entails absorbing a 
quantum of flux (k = 1) while following I implies no change in the flux trapped by the JJN. 
Of course, in equilibrium the system will prefer to follow 77 rather than /, which is a higher 
energy state. However, if there is an external f e {r) present the JJN will be allowed to follow 
curve / as well. We shall call the process following curve II a phase slip (PS) case while the 
one that follows I the no phase slip (NPS) process. For the £ ^ case, the ground state 
energy becomes = —cos[^-([f]k — /)] + and the corresponding normalized current is 
%n = sin[^([/]fc — /)]. From these equations one finds that the total flux in the ring is given 
by / = fe + 4a s ^ n [lf ([/]* ~~ /)]' which is a self-consistent transcendental equation that, for 
a given set of values for /° and £, can be solved numerically. We note that in contrast to 
the N — 1 case, where there is a critical value £ c — l that separates non-hysteretic (I > 1) 
from hysteretic (I < 1) behavior, the JJN with iV > 2 is always hysteretic. The size of 
the hysteresis loop grows as £ increases from zero until £ > N, when the hysteresis loop 
covers the whole /°e[0, 1] range. The metastability of the JJN is evident from the history- 
dependent ground state energy shown as a dashed line in Fig. 1. In the PS case we see 
that there are two values of /° for which has a discontinuity, one while ramping up 
the other when decreasing it. The specific / e ° values at which the discontinuities take place 
depend on £ and N. The NPS case is, again, represented by continuous processes. 

With the above information we are now ready to discuss the dynamic properties of 
the JJN. We start by introducing the important parameter k 2 = v^jv^ with v<& the 
characteristic frequency for the relaxation of magnetic fields. In the linearized regime of the 
Josephson term in Eq(3), v<$, = R/L, and k = 2 ni l ^ \A/^' P arame ^ er measures 

the importance of screening effects: when k — oo they are negligible while for k < 1 they 
are all important. Thus we want to study the dynamics of the JJN as a function of k in the 
two extreme regimes of PS and NPS. For simplicity we shall write £ instead of k in our 
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results, although the connection between k = yl/l is only valid in the linearized regime. 
Most of the time the value of [/(r)]fc is constant except when /(r) passes through the phase 
slip points of Fig. 1. In this case [f(r)]k is discontinuous and its derivative -^[f{r)]k has 
a 5-function character. We use this information in evaluating the time-averaged current 
given in Eq(3). We start by considering the simplest case when I = 0, i.e. when there is no 
screening. In the NPS case < ijq > can be calculated analytically since ^"[/(t)]^ = 0, and 
the integral < i N >= lim^^ \ f__ 2 r ds sm[^(2irk - 27r/ e ° + v dc s + sin 27n^s)] + v dc /N, 
gives < i N >= v dc /N, when v dc ^ mNu^, and < i N >= sin[^(A; - fe)) J -m(^^-) + v dc /N, 
whenever v dc = mNv^, with m — 0, 1, 2, • • -. In deriving this result we have used the identity 
gj^sm^ _ Yf^oo J -m{z)e~ vm ' t> , with J - m {z) the Bessel function of integer order. Note that 
when v dc = mNu^ the initial condition /° plays a crucial role since < zjv >=< «jv > (/e); 
and thus when frf c = mNv^ there are finite intervals of < > for each value of v dc . 
These giant Shapiro resonances (GSR) are shown in Fig. 2(a) for iV = 10, = 0.01 and 
v ac = 1- We note that the size of the current intervals for which there are resonances are 
comparable to one. The averaged current in the PS can be evaluated analytically following 
a similar logic as in ref . . The qualitative result is that there are subharmonic resonances 
for v dc = ^Nu^, with n = 1,2 ■ ■ ■. The analytic analysis does not give, however, the order of 
magnitude of these subharmonic resonances. To resolve this question we evaluated < in > 
in the PS case numerically. The results for iV = 10, v$ = .01, v ac = 1 are shown in Fig. 2(c). 
These results were obtained as follows: first the time interval was divided in At segments, 
with i/^Ar = 10~ 4 , and the v dc interval into 24 pieces. Next, for each value of v dc we took 100 
initial conditions for /° . Whenever / e (r) reaches a discontinuity point, which is determined 
with a precision of 10~ 6 , we changed the value of [f(r)}k- Note that the scale in Fig. 2(c) is 
about an order of magnitude smaller than in Fig. 2(a). 

We now consider the general case with i ^ 0. To evaluate the time-averaged current 
in this case we need to solve the differential equation given in Eq(3). This is not a simple 
ordinary differential equation since there is an implicit transcendental dependence of ijv(r) 
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through the sine function. The equation is solved by using a fourth order Runge-Kutta 
method following a similar procedure as in the £ = case, for parameter values similar 
to those of Fig. 2(c) but with £ = 50.0, or k — 0.1414. For the NPS case the results 
are shown in Fig. 2(b). There we see subharmonic resonances but of magnitude larger 
than in Fig. 2(c). We call these resonances then giant subharmonic Shapiro resonances 
(GSSR). In the PS case and for the same parameter values as in Fig. 2(b), Fig. 2(d) 
shows, in contrast, that the subharmonic resonances remain of the same order of magnitude 
as in the £ = limit. We also calculated the spectral function, not shown here, defined 
as S(27iu l j ) ) = Hindoo \^Jq i^^e^^als] 2 . It is found that at a resonance the current as a 
function of time is a periodic function but with a complicated subharmonic structure. 

The results presented in Fig. 2 show that the JJN is capable of exhibiting truly GSSR 
in the NPS case although they are also present, but of a smaller magnitude, in the PS 
case. To further understand the properties of these resonances we calculated the magnitude 
of the resonance spikes widths A < i N > as a function of \fi = 1/K,v ac and N. The 
corresponding results are shown in Fig. 3(a-d) for the 1/2 resonance, for both the NPS 
case (o) and the PS case (x). The significant characteristic of Figs. 3(a-c) is that for the 
NPS case there are maximum or optimal A < ijy > values as a function u^, VI and N. The 
situation is less clear for the PS case although small maxima can be seen in Figs. 3(a-c). 
These specific maximum values are listed in the figure caption. Note that as a function of 
£, and as seen in Figs. 2, A < i N > is larger for the NPS case in the large £ limit but it 
is slightly smaller in the PS case in the small £ regime. A similar suppression of the 1/2 
step width for large was found in [3,7] for a current driven square array. The suppression 
was not found, however, for integer steps. In contrast, we found that the n — 1 step width 
as a function of behaves in a qualitatively similar way as in the 1/2 case. In Fig. 3(d) 
we show the behavior of A < i N > as a function of v ac . Contrary to what is seen in one 
J J or in the 2D-array calculations, A < i^ > grows monotonically up to about v ac = 3 and 
then it appears to become periodic for the NPS case but is still aperiodic for the PS limit. 
We checked that these results are stable against an increase in the number of periods used 
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to calculate the averages. All the results shown in Fig. 3 indicate that the qualitative and 
quantitative properties of the JJN are inherently different from those obtained previously 
for one JJ or for 2D-arrays. 

These results allow us to make a rough estimate of the appropriate experimental param- 
eters for which these resonances can be seen experimentally taking the typical values from 
experiments HHI. For a lattice constant a ~ lOfim and with separation between junctions 
of b ~ 2/^m, with the optimal value of N = 16, we get a diameter of d ~ 61/zm. The critical 
currents are between 1 x 1CT 2 to 10 mA, so that using the estimate for the inductance for 
a ring |J of L = 1.25/iod, one gets the values of \/J between 1.73 ~ 54.84 which are within 
the range of maximum values shown in Fig. 3(b). 

In conclusion, we have presented a detailed analysis of the response of a Josephson 
junction necklace to an external time-dependent flux. The main result of our analysis is 
that there are giant subharmonic Shapiro resonances, mainly in the no phase slip or "fast" 
regime. Furthermore, we found that these resonances have a strong dependence on the 
number of junctions N, normalized frequency v$ and effective inductance i. An estimate 
was also provided for the possible experimental conditions under which these resonances 
may be seen. 
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FIGURES 

FIG. 1. Normalized Ground state energy for N = 10 as a function of /° . Solid line £ = and 
dashed line £ = 5.0. The curves I and II represent the no-phase-slip (NPS) and phase-slip (PS) 
processes, respectively. The inset gives a schematic representation of the JJN model studied here. 

FIG. 2. The time-averaged normalized current < in > vs Vdc for N = 10, v§ = 0.01 and 

v ac = 1-0. (a) the iVPS and (c) PS cases, both with £ = 0. (b) the NPS and (d) PS cases with 

£ = 50.0. The results were obtained as described in the text. The — and — subharmonic resonances 

2 3 

are indicated by arrows. 

FIG. 3. Renormalized current width A < ijy > for the 1/2-resonance as a function of (a) 
z/0, (b) y/I, (c) N and (d) v ac in the NPS (o) and PS (x) cases. In (a) v ac = 1,7V = 10 and 
£ = 50.0. In (b) i/0 = 0.01, N = 10 and w ac = 1. In (c) w ac = 1,4 = 50.0 and = .01. In 
(d) i/0 = 0.01, £ = 50.0 and N = 10. The approximate maxima in Figs. 3(a-d) are located at 
u rnax ^ q.0251, yj l max ~ 5.012 and N max ~ 16, respectively. 
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